#include<iostream>

using namespace std;
typedef long long ll;
const int N=1e7+10;
bool st[N];
int p[N],cnt,phi[N];
ll f[N];
void get_phi(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!st[i])
		{
			p[cnt++]=i;
			phi[i]=i-1;
		}

		for(int j=0;1ll*i*p[j]<=n;j++)
		{
			int x=i*p[j];
			st[x]=true;

			if(i%p[j]==0)
			{
				phi[x]=phi[i]*p[j];
				break;
			}
			else
			{
				phi[x]=phi[i]*(p[j]-1);
			}
		}
	}
}
int main()
{
	int n;cin>>n;
	get_phi(n);
	for(int i=1;i<=n;i++)
	{
		f[i]=f[i-1]+phi[i];
	}

	ll ret=0;
	for(int i=0;i<cnt;i++)
	{
		ret+=f[n/p[i]]*2+1;
	}

	cout<<ret<<endl;
}